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ABSTRACT 

We report the discovery of burst oscillations from the intermittent accretion-powered millisecond 
pulsar (AMP) HETE J1900. 1-2455, with a frequency ^ 1 Hz below the known spin frequency. The 
burst oscillation properties are far more similar to those of the non-AMPs and Aql X-1 (an intermittent 
AMP with a far lower duty cycle), than those of the AMPs SAX J1808.4-3658 and XTE J1814-338. 
We discuss the implications for models of the burst oscillation and intcrmittency mechanisms. 
Subject headings: binaries: general, stars: individual (HETE J1900. 1-2455), stars: neutron, stars: 
rotation, X-rays: bursts. X-rays: stars 
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1. INTRODUCTION 

There are two ways in which we can measure the spin of 
rapidly rotating accreting neutron stars: via accretion- 
powered pulsations (due to magnetic channeling of in- 
falling material), or burst oscillations (coherent pulsa- 
tions seen during the thermonuclear explosions of ac- 
creted fuel that give rise to Type I X-ray bursts). The 
mechanism that generates the brightness asymmetries re- 
sponsible for burst oscillations is not yet known. Once 
understood, however, it should shed light on the ther- 
monuclear burning process and the composition of the 
outer layers of the neutron star (jStrohmaver fc BildstenI 
[200llWatti[200l . 

Burst oscillations often drift upwards in frequency 
by a few Hz during a burst. However the stabil- 
ity of the asymptot ic frequency for a particular source 
uno et all l2002af ). and the detection of the same 
frequencies in superbursts as w ell as Type I bursts 
(jStrohmaver fc Markwardt]l2002t l. suggested even at an 
early stage that burst oscillation frequency might be 
(to within a few Hz) the spin frequency. This was 
confirmed in 2003 with the discovery of burst oscilla- 
tions at or within a few Hz of the spin frequency in 
the accretion-powered millisecond pu lsars (AMPs) SAX 
J1808.4-3658 and XTE J18 14-338 (jChakrabartv et all 
[2003tlStrohmaver et al.l[2003h . The near identity of burst 
oscillation frequency and spin frequency is now the major 
constraint on models for the burst oscillation mechanism. 

The AMPs are particularly useful when investigating 
the burst oscillation mechanism, since they are the only 
sources in which we can gauge the effects of the mag- 
netic field and asymmetric fuel deposition. They are 
also the only stars where we can measure the size and 
sign of the small offset between the spin and burst oscil- 
lation frequency. However both SAX J1808. 4-3658 and 
XTE J1814-338 have burst oscillations with quite atyp- 
ical properties compared to the other sources (Section 
[3]). It is not yet clear whether we are seeing a contin- 
uum of behavior (which could be explained by one burst 
oscillation mechanism) or completely separate classes. 

The picture has evolved in the last three years 
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with the discovery of three intermittent AMPs: 
HETE J1900 1-245 5 SAX J1748.9-2021 and Aql X-1 
(iKaaret et al.ll2006l: iGalloway et ani2007l: iGavriil et all 
I2007t I Altamirano et al.ll2008l : ICasella et al.ll2008[ ). Un- 
like the other AMPs these sources show accretion- 
powered pulsations only sporadically. What cause the 
intcrmittency is also not known, with mod els including 
spora dic obscuration of the magnetic poles (|Gogiis et al.1 
l2007l and references therein) , wandering of the accretion 
funnel and hence the hot spot around the mag netic pole 
(|Romanova et al.ll2003l.l2004l :'Lamb e t al.ll2008h. or mag - 
netic field burial (jCumming e t al. 20011: ICumminel[2008f ). 

If the magnetic field is what makes the burst oscilla- 
tions of the AMPs so different from those of the other 
stars, then an intcrmittency mechanism that affects the 
magnetic field may also influence burst oscillation prop- 
erties. Burst oscillations from the intermittent AMPs 
(all of which burst) are therefore an important piece of 
the puzzle. Until now, however, only Aql X-1 has shown 
burst oscillations, at ~ 0.5 Hz below t he spin frequency 
(jZhang et al.lll998t ICasella et al.ll2008f) . The properties 
of the Aql X-1 burst oscillations are very similar to those 
of the non-AMPs: however this source has the lowest 
pulsation duty cycle of any of the intermittent AMPs. 
In this Letter we report the discovery of burst oscilla- 
tions from the int ermittent AMP HETE J1900. 1-2455 
(jWatts et al.ll2009l ) , and consider the implications for the 
burst oscillation and intcrmittency mechanisms. 

2. BURST ANALYSIS 

HETE J1900. 1-2455 was first detected in June 2005, 
and has remained in outburst eve r since apart from a 3 
week period of quiescence in 2007 (jDegenaar et al.|[2007f) . 
It was quickly identified as an AMP with a spin frequency 
of 377.3 Hz, an orbital period of 83.3 minutes, an d a low 
(< O.IMq) mass companion (jKaaret et al.ll2006l ). It is 
at a distance of « 5 kpc, determined from Eddington 
limited X-ray bursts (jCallowav et al.ir2008af l. Accretion- 
powered pulsations of low fractional amplitude (< 3 %) 
were detected sporadically during the firs t two months of 
outburst, bu t have not been seen since (jCallowav et al.1 
120071 l2008bf ). There is some evidence that the in- 
termittent appearance of the accretion-powered pulsa- 
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tions may be related to the occurrence of Type I X-ray 
bursts, but whether th e relationship is causal is not clear 
()Gallowav et al.ll2007[ ). From 2005-2008 several X-ray 
bursts from this source were detected by HETE-II, the 
Rossi X-ray Timing Explorer (RXTE) and SWIFT^. We 
searched the bursts detected by both the RXTE Propor- 
tional Counter Array (PCA) and the SWIFT Burst Alert 
Telescope (BAT) for oscillations, but found no significant 
signal. 

2.1. Burst oscillations 

During routine monitoring on April 2 2009 (08:57 
UTC), RXTE detected another X-ray burst with both 
the PCA and the High Energy X-ray Timing Experiment 
(HEXTE). Timing analysis was conducted using 125 /iS 
time resolution PCA event mode data from the two active 
PCUs (PCUs and 2). We used aU photons in the 2-30 
keV range, the band where burst emission exceeds the 
persistent level. The data were barycentered using the 
JPL DE405 and spacecraft ephemerides, with the source 
position of Fox (2G0_§). Some event mode data overruns 
occurred in the burst peak, leading to short data gaps. 

A dynamical power spectrum (Figure [1]) reveals strong 
burst oscillations during the initial decay of the burst, 
drifting upwards by about 1 Hz. Selecting only data af- 
ter t=4 s in Figure [1] (the time of the final data gap), 
wc find maximum Leahy power in the range 17-30 in 4 
independent consecutive 2 s bins, an extremely robust 
detection. The oscillations fall below the detectability 
threshold when the frequency is 376.3 Hz, 1 Hz below 
the spin frequency. Although the binary ephemcris can- 
not be extended to 2009 (the errors are too large), orbital 
Doppler effects will shift th e spin frequency by 0.009 Hz 
at most (jKaaret et al.ll2006[ ). so the offset in the frequen- 
cies is secure. A search for accretion-powered pulsations 
during the observation revealed no significant signals. 

To study the pulse profile, we extracted data for 
the 8s period (t = 4- 12s in Figure [J) during which 
the oscillations are detected. We fitted the frequency 
drift using a polynomial model to maximise the power. 
Power was maximised with a quadratic drift model, 
V = vq i>Q{t ~ to) + h'o{t — Iq)'^ , which yielded a Leahy 
normalised power of 88 (the best linear model gave a 
maximum power of 78, and higher order polynomial 
terms yielded no further improvement). Using the best 
fit quadratic frequency model we generated a folded 
profile. The persistent (non-burst) emission level was 
estimated from the ~ 1000s prior to the burst. The 
resulting pulse profile is shown inset in Figure [TJ The 
pulse is well fit with a simple sinusoid (x^/ dof~ 31/37), 
with no need for harmonic content. The fractional 
amplitude was (3.5 ± 0.3)% Root Mean Square (RMS) 
(2-30 keV), comparable to the fractional amplitude of 
the accretion-powered pulsations of this source. The 
amplitude of the pulsations rises with energy, from 
(3.1 ± 0.5) % RMS in the 2-5 keV band to (5.5 ± 0.6)% 
RMS in the 10-20 keV band. The burst oscillations also 
show hard lags, with the 10-20 keV pulse lagging the 2-5 
kcV pulse by (0.09 ± 0.03) cycles. 

^ IGallowav et al.l l)2008bl ') list most of the bursts detected during 
2005-2008. One additional burst during this period was recorded by 
RXTE on December 5 2007 (05:57 UTC), and there is a candidate 
burst in slew data on November 27 2007. 
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Fig. 1. — Lightcurve and dynamical power spectrum (2-30 keV) 
for the April 2 2009 burst. The dynamical power spectrum uses 
overlapping 4 s bins, with new bins starting at 0.25 s intervals. 
We use a Nyquist frequency of 2048 Hz and an interbin response 
function to reduce artificia l drops in a mplitude as the frequency 
drifts between Fourier bins l lvan der Klisll989f l . The contours show 
Leahy normalized powers of 20-40, increasing in steps of 5. The 
dashed line indicates the 377.3 Hz spin frequency. The inset shows 
the pulse profile of the burst oscillations for the 8s interval whore 
they are detected, folded with the best fit polynomial frequency 
model (counts per 0.05 cycle phase bin against phase in cycles), 
with the mean subtracted. The resulting profile is well fit by a 
simple sinusoid (overplotted) with no requirement for any harmonic 
content. Two cycles are shown for clarity. 

2.2. Burst properties 

To study the characteristics of the burst containing the 
oscillations, we extracted spectra every 0.25s from the 
PCA Event data (E_125us_64M_0_ls), which has 64 en- 
ergy channels between 2 and 60 keV. We used data from 
PCUs and 2, which were operating at the time of the 
X-ray burst. We generated an instrument response ma- 
trix for each spectrum, and fitted it using XSPEC version 
11.3.2. We added a 1% systematic error to the spectra 
and restricted the spectral fits to the energy range 2-25 
keV. We extracted spectra using 100s segments before 
and after the burst, and used them as background in 
our fits. We fi tted each of the 0.25s spectra with an ab- 
sorbed (wabs. [Morrison fc McCammonlll983[ ) blackbody 
model (bbodyrad model in XSPEC). This approach, 
which takes account of the contribution of the persistent 
and background emission, is sta ndard procedure in X -ray 
burst analysis (see for example iKuulkers et al.|[2003f ). 

The simple blackbody fit described above yields a peak 
luminosity of 3.6 x 10'^* ergs s~^, and a total energy 
of 2.7 X 10'^^ ergs. The evolution of the temperature 
Thh and the radius Rbb given by the spectral fits resem- 
bles th ose seen in photospher ic radius expansion (PRE) 
bursts (jCallowav et al.|[2008aD . i.e. (i) Rbb reaches a local 
maximum close to the time of the peak flux, (ii) lower 
values of Rbb were measured following the flux maximum 
and (iii) there is a local minimum in Tbb at the same time 
as the maximum in Rbb- This conclusion (and the de- 
rived peak flux and fluence) should be treated with care, 
however, since the x^/dof of the spectral fits can reach 
values as high as 4 or 5. The problem stems from an 
excess at high energies . Similar issues were reported by 
IGallowav et al.l (|2008b[ ) for other bursts from this source, 
and the physical cause is not yet understood. Adding a 
power law to our model, for example, reduces x^ /dof to 
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Fig. 2. — 2-25 kcV lightcurves for the bursts recorded by the 
RXTE PCA from 2005-2009. Background and persistent emission 
have been subtracted and deadtime corrections applied (this was 
not done in Figure [TJ. The bursts are numbered according to their 
position in the color-color diagram (Figure |4]l. B6, the burst with 
the oscillations, has quite a different lightcurve. 

below 2: the evolution of Ti,b and Rbb in this case still 
suggests a PRE burst. Photospheric touchdown occurs 
3.5 — 5s after the start of the burst, consistent with the 
first detection of the burst oscillations. 

The parameters of the April 2 burst do not differ dra- 
matically from those of other bursts from HETE J 1900, 
which also show evidence for PRE subject to unc er- 
tainties about the spectral fits (jGallowav et al.|[2008al lbl) . 
The burst with oscillations does not however have the 
extended or double peak structure exhibited by most of 
the other bursts (Figure [2). Th e burst lightcur v e show s 
a fast rise, « 0.4 s defined as in iGallowav etldl (|2008a[ ). 
The decay can be modelled with a double exponential 
with decay timescales as of 7.3 s and 8.4 s respectively. 
Total burst duration is w 60 s (burst end is defined as 
the time when flux falls to 1% of the peak flux, corrected 
for persistent emission). 

2.3. Persistent emission 

Figure [3] shows the intensity, defined as the coun- 
trate in the 2.0-16.0 keV band, for all publicly avail- 
able observations from 2005-2009. Figure [J] shows the 
corresponding color-color diagram. Data were deadtime 
corrected, background subtracted, and X-ray bursts re- 
moved. For each observation we calculated X-ray colors 
from the Standard2 data. We defined soft color as the 
ratio between count rates in the 3.5-6.0 and 2.0-3.5 keV 
bands and hard color as the ratio between count rates 
in the 9.7-16.0 and 6.0-9.7 kcV bands. We normalized 
colors and intensity t o the Crab values nearest in time 
(jKuulkers et al.lll994D and in the same P CA gain epoch 
(see for example lvan Straaten et aLll2003[ ). 

It is clear from Figures that the observation where 
the burst oscillations were detected was rather unusual. 
The burst happened when source intensity was at its 
highest recorded level (w 63 mCrab) and when it was 
in the soft (banana) state. All previous bursts have been 
detected in harder states. Fitting the PCA spectrum 
with an absorbed disk-blackbody plus power-law model, 
a nd assuming a standar d bolomctric correction factor of 
/ flm 't Zand et al]|2007[ ) we find an unabsorbed bolomet- 
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Fig. 3. — Average 2-16 keV intensity per RXTE observation for 
HETE J1900. 1-2455 from 2005 to April 4 2009 (data from March- 
October 2008 are not yet public). Observations containing bursts 
are marked in black: triangles (no oscillations), square (oscilla- 
tions). We include the four confirmed bursts detected by the PCA 
from 2005-2008 (the candidate slew burst also occurs in the hard 
state, but is not shown), the burst with oscillations from April 2 
2009 and an additional burst detected by the PCA on April 4 2009. 
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Fig. 4.— Color-color diagram for HETE J1900.1-2455 using aU 
publicly available RXTE data from 2005 to April 4 2009 (excluding 
observations with intensities below 2 m Crab, see Figure [3ll. The 
sourc e behaves like a typical atoll source l IHasinger &: van der Klii 
119891) . The grey dots show the colors computed from 160s segments 
of data. Observations containing bursts are marked in black: tri- 



angles (no oscillations), square (oscillations), 
bursts are shown in Figure [2] 
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tion was fixed to 1.6 x 10 cm ~ ). At a distance of 5 
kpc, this corresponds to 4% of the Eddington luminosity 
if wc assume iEdd = 2.5 x 10^^^ erg s~^. 

3. DISCUSSION 

What causes burst oscillations is still not understood. 
Flame spread from a point should lead to asymmetries 
in the early phase of the burst. However, while this may 
explain the presence of burst oscillations during the rise, 
it is not sufRcient to explain the continuing presence of 
large-scale asymmetries in the burst tail once the en- 
tire s tellar surface has ignited (jStrohmaver fc BildstenI 
|2006[ ). This has led to the consideration of alternative 
models. Non-radial global oscillations in the surface lay- 
ers of the star, excit ed by the flame spread, remain a 
promising possibility ()Hevlll200^ . In this case the bright- 



ric flux of 3.5 x 10 ergs s cm (interstellar absorp- ncss asymmetry would be caused by variations in ocean 
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height associated with the mode. Inertial frame pattern 
speed can be very close to the spin rate, and coohng 
of the layers in the aftermath of the burst would nat- 
urally lead to frequency drift. The major problem of 
this model is that it overpredicts the size of the fre- 
quen c y drift compared to obse rvations (jPiro fc BildstenI 
120051: iBerkhout fc Levinll2008|) Alt ernative possibilities 
like ph otospheric modes ( Hevlll200^ or shear oscillations 
()Cumm ing 20 0^ also have sho rtcomings in their present 
form ( Berkhout fc Levinir2008f ). Current efforts to resolve 
these problems are focusing on the role of the magnetic 
field, which can dominate the dynamics in the surface 
layers and is expected to have a large effect on surface 
modes. 

If the magnetic field is the most important factor deter- 
mining the properties of burst oscillations then this might 
lead to a natural explanation for the differences in the 
properties of the burst oscillations of SAX J1808. 4-3658 
and X TE J1814-338 (hereafter collectively referred to 
as SX: ICha krabartv ct al."2003': 'Strohmavcr ctalj[200l; 
IWatts et al.. 2005 : Watts fc Strohmayer 2006) c om pared 



a- /yub,] c om pared 
ai.l l2002allbl . ^003l . 



to those of the non-AMPs jMuno et 

120041 iGallowav et"an l2008af ) . These differences can be 
summarized as follows. (1) SX show oscillations in all 
of their bursts, even in the hard state; the non-AMPs 
do not (oscillations are detected primarily although not 
exclusively in the soft state). (2) SX show oscillations 
throughout their bursts (except during PRE); the non- 
AMPs, in most cases, do not. (3) SX burst oscillations 
have amplitudes that fall with energy; the non-AMP os- 
cillations have amplitudes that rise. (4) SX burst oscilla- 
tions have detectable harmonic content; the non-AMPs 
have none. (5) Frequency drift in SX is either very fast 
or non-existent; in the non-AMPs it is much slower. 

The magnetic field also plays an important role in mod- 
els for intermittency of accretion-powered pulsations. In 
the obscuration and accretion stream wander models, the 
field is always present at the level necessary to channel 
the flow: but the accretion hot spot either wanders out 
of the line of sight or is obscured by magnctosphcric ma- 
terial. In the magnetic burial model the field strength 
and geometry change as the field is suppressed by ac- 



cretion. If magnetic field affects both intermittency and 
burst oscillations, then studying the latter may enable us 
to pinpoint the cause of the former. 

In Aql X-1 the burst oscillations are similar in proper- 
ties to those of the non-AMPs: however this source has 
an exceptionally low duty cycle so its intermittent pul- 
sation episode may have been triggered by an extremely 
rare event. In HETE J1900. 1-2455 the accretion-powered 
pulsations lasted for much longer. Its burst oscillations, 
however, also behave like those of the non-AMPs. This 
has important implications. If hot spot wander or ob- 
scuration models for intermittency are correct, then a 
field strong enough to channel infalling material cannot 
be the only factor causing the atypical burst oscillations 
of SAX J1808.4-3658 and XTE J1814-338. Some other 
factor, such as the presence of a strong temperature gra- 
dient around the magnetic pole, niust also play a ro le in 
the burst oscillation mechanism ijWatts et al.ll2008[ ). If 
on the other hand the field has been buried in HETE 
J 1900. 1-2455, then it must be screened to a depth where 
it is unable to affect the burst oscillation mechanism on 
the timescale of the burst. Detailed calculations will be 
required to resolve this issue, but if screening at the burn- 
ing depth is not possible, then this could point to pho- 
tospheric rather than oceanic modes as a cause of burst 
oscillations. 

One other point of note is that in both Aql X-1 and 
HETE J1900. 1-2455, burst oscillation frequency remains 
below spin frequency. Given the similarities in properties 
it seems probable that this is the case for non-AMPs as 
well. However it is not the case for SAX J1808.4-3658 
and XTE J1814-338. In SAX J1808.4-3658 the burst os- 
cillations are first detected below the spin frequency but 
rapidly overshoot it, settli ng ~ 0.1 Hz above the spin 
frequency in the burst tail (jChakrabartv et ani2003f ). In 
XTE J1814-338 the two frequencies are identical, and in 
fact burst oscillations and accre tion-powered pulsations 
are coherent and phase -locked (jStrohmaver et all l2003l : 
IWatts et al.1l2005ll2008[ ). Any unified model of burst os- 
cillations must be able to explain this diversity in the re- 
lationship between burst oscillation frequency and spin. 
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